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Abstract 

Wc report the first results of a near- and mid- infrared deep survey with the Infrared Camera (IRC) 
onboard AKARI in the performance verification phase. Simultaneous observations by the NIR, MIR-S and 
MIR-L channels of the IRC with effective integration times of 4529, 4908, and 4417 seconds at 3, 7, and 
15 fiui, covering 86.0, 70.3, and 77.3 arcmin^ area, detected 955, 298 and 277 sources, respectively. The 
5 a detection limits of the survey are 6.0, 31.5 and 71.2/xJy and the 50 % completeness limit are 24.0, 
47.5, and 88.1/xJy at 3, 7, and 15 fim, respectively. The observation is limited by source confusion at 3 
jum. We have confirmed the turnover in the 15 /um differential source counts around 400 fiJy, previously 
detected by surveys with the Infrared Space Observatory. The faint end of 15 fim raw source counts agree 
with the results from the deep surveys in the GOODS fields carried out with the Spitzer IRS peak up 
imager and the predictions of current galaxy evolution models. These results indicate that deep surveys 
with comprehensive wavelength coverage at mid-infrared wavelength are very important to investigate the 
evolution of infrared galaxies at high redshifts. 
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1. Introduction 

The Infrared Camera (IRC; Onaka et al. 2007; Wada 
et al. 2002) is a near- and mid- infrared scientific instru- 
ment onboard the first dedicated Japanese infrared astro- 
nomical satellite AKARI (Murakami 2007). The unique 
instrumental features of the IRC, compared to the simi- 
lar contemporary IRAC instrument on the Spitzer Space 
Telescope (Fazio et al. 2004), are a wider field of view 
(FOV; 10 X 10 arcmin^) and a more comprehensive con- 
tiguous wavelength coverage (2-26/im) with a larger num- 
bers of filters (9 bands) of comparable sensitivity to IRAC. 
A deep and wide field near- and mid- infrared survey 
(LS NEP DEEP) covering these wavelengths by all the 
9 bands is underway in the North Ecliptic Pole (NEP) 
region (Matsuhara et al. 2006). 

Previous deep surveys at mid-infrared wavelengths have 
made the striking discovery of a bump in the Euclidean- 
normalized differential source counts around 0.4 mJy at 15 
fim in the ISO/ISOCAM surveys (Elbaz et al. 1999) and 
0.3 mJy at 24 fim in the Spitzer MIPS surveys (Papovich 
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et al. 2004). Mar lean et al. (2004) compared the 15 ^m 
and 24 fj. source counts adopting a 24/15 fim flux ratio of 
1.2, which was derived from the SED template of Chary 
& Elbaz (2001). They found that a shift in the turnover 
of the Euclidean-normalized counts, however, the fainter 
side of 15 fim counts was not well constrained due to the 
limitation of the sensitivity and areal coverage of the ISO 
surveys. Pearson (2005) showed that models incorporat- 
ing either a burst of evolution or a more smoother contin- 
uous evolution could simultaneously reproduce both the 
15 /im ISO and 24 /im Spitzer counts and concluded that 
revealing the dominant population at the fainter end of 
the 15 fim counts (e.g., starburst, LIRGs, ULIRGs) was 
an important point to be resolved. 

One of the objectives of the LS NEP DEEP survey is to 
reveal the nature of the faint 15 fim population through 
a deep and wide multi-band survey together with the 
SED fitting technique focused on the PAH emission and 
the silicate absorption features (Takagi & Pearson 2005). 
Note that although Spitzer also has the capability, via its 
IRS 16fim peak-up Imaging (PUI) mode, for imaging in 
this wavelength range (Teplitz et al. 2005a; Kasliwal et 
al. 2005), however, the FOV of the IGfim PUI mode is 
very small (1 x 1 ; 100 times smaller compared to that 
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of the IRC), and moreover, there still exists a gap in the 
wavelength coverage between the IRS 16/im PUI and the 
shorter 8/im IRAC band. 

In order to evaluate the capability of such a deep and 
wide survey by the IRC, we have conducted a pilot sur- 
vey in both the imaging and spectroscopic modes with 
10 pointed observations for each mode at the end of the 
performance verification phase of the AKARI satellite. 

In this paper, we report on the first results of the pi- 
lot imaging survey. The result of the corresponding pi- 
lot spectroscopy survey will be reported elsewhere. In 
this work, section 2 describes the observations. In section 
3, the data reduction is presented, source extraction and 
photometry are described in section 4, section 5 describes 
the estimation of sensitivity and completeness of the sur- 
vey. Finally, we discuss the capability of the IRC in deep 
survey mode and the 15 fim source counts in section 6. 

2. Observation 

2.1. AOT for deep surveys 

Since the main objective of the AKARI mission is the 
far-infrared All-Sky Survey (Pearson et al. 2007), the or- 
bit is by design a Sun-synchronous polar orbit which puts 
severe constraints on the visibility of individual target 
fields on the sky. Thus, the operation for any pointed 
observations have been made as simple as possible, pro- 
ducing strong constrains on the Astronomical Observation 
Templates (AOT) for the IRC pointed observations such 
as (1) fixed frame integration time, (2) no mapping mode 
in a single pointed observation (approximately 10 min- 
utes), that were not apparent for infrared space observa- 
tories such as Spitzer or the Infrared Space Observatory 
(ISO)i. 

A series of three AOT sets for pointed observations with 
the IRC have been prepared. The first set, referred to as 
IRC02 and IRC03, are for medium deep multi-band sur- 
veys, in which both changes of filters and dithering in 
target position is performed. The second set are for deep 
surveys and are referred to as IRCOO and IRC05, in which 
no filter change and no dithering is performed, in order 
to minimize the dead time and maximize the observation 
time. For the IRCOO and IRC05 AOTs, at least three 
pointed observations are required in order to ensure re- 
dundant and reliable observations. The third AOT set, 
referred to as IRC04 is designed for spectroscopic obser- 
vations and is described elsewhere (Ohyama et al. 2007). 

Using the AOTs designed for deep surveys, a deep 
and wide field near- and mid- infrared survey (LS NEP 
DEEP) is underway in the North Ecliptic Pole (NEP) re- 
gion (Matsuhara et al. 2006). In order to evaluate the 
capability of deep imaging surveys by the IRC, we have 
conducted a pilot survey with 10 pointed observations. 

^ These constraints make the IRC inefficient for wide stiallow sur- 
veys in which a depth of one pointed observation or better is 
not required. We have prepared a separate IRC scanning mode 
for such a shallow and wide surveys including the IRC All-Sky 
Survey at mid-infrared wavelengths (Ishihara et al. 2006). 



2.2. pilot survey 

The observations were executed from April 29 to May 
6 in 2006, almost at the end of the performance verifi- 
cation (PV) phase of AKARI. We selected a target po- 
sition of (alpha, delta)=(268.88, 66.61), chosen by the 
following reasons. (1) The FOV of the NIR and MIR- 
S channels must be within the continuous viewing zone of 
AKARI (CVZ: within 0.6 degree from the North or South 
ecliptic poles). (2) The field must be well studied by the 
other observations in order to make cross-identification 
of the AKARI sources. The deep multi-color optical im- 
ages (Wada et al. 2007) taken by Subaru/Suprime-Cam 
(Miyazaki et al. 2002) were used for this purpose. (3) The 
field must have few bright objects in order to make the 
deep survey. We carefully selected a field, which had no 
bright stars, in the Suprime-Cam image field. 

We used the AOT for deep observation (IRCOO) with 
an AOT filter combination parameter of "b" , which cor- 
responds to the filter set of the N3 (3.3 /xm), S7 (7.0 /xm) 
and L15 (15.0 yum) band filters, for the NIR, MIR-S and 
MIR-L channels, respectively. Wc used a target position 
parameter of "N" , which means that the center of the FOV 
of the NIR channel (and simultaneously, the MIR-S chan- 
nel which observes the same point on the sky via a beam 
splitter) targeted the specified coordinates, while that of 
the MIR-L channel was pointed at a parallel field, around 
25 arcmin apart from that of the NIR/MIR-S. No dither- 
ing operation between pointed observations were made, 
i.e. we used the same target position for all pointed ob- 
servations^. However, the FOV of the NIR/MIR-S rotates 
around the center of the field as the direction of the Sun 
changes (approximately one degree per day), resulting in 
good redundancy except at the very center of the FOV. 
Note that good redundancy in the observations is possi- 
ble over the entire FOV of the MIR-L channel, since the 
FOV of the MIR-L rotates around the FOV of NIR/MIR- 
S. Details of the observation parameters are summarized 
in Table 1. 

Typically, 10 and 30 frames with effective exposure 
times of 44.4 and 16.4 seconds in each frame were ob- 
tained for a single pointed observation, for the NIR and 
MIR-S/L channels, respectively. The total exposure time 
was 4530, and 4908 and 4417 seconds, for the NIR, MIR-S 
and MIR-L channels, respectively. 

3. Data reduction 

The data were reduced by the standard IRC imaging 
pipeline (version 060801; see IRC Data User Manual, 
Lorente et al. 2007). The dark frame (version 060428), 
which was measured just before the PV phase observations 
began, was subtracted. Flat fielding was performed based 
on the flat field frame (version 060626), which was created 
from 6 pointed observations of a high background region 
near the ecliptic plane with the AOT IRC02. Distortion 
correction (version 060529) was performed based on the 

^ Dithering operation between pointed observations have already 
implemented in the observations after 2006-09-24X09:00:00 
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Table 1. Observation log 
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Target id 


AOT 




OBb 


-DATE 




N/S position 


L 


position 
















RA 


DEC 


RA 




DEC 


5020053.1 


IRCOO 


b;N 


2006- 


05- 


03103 


32:41 


268.88 


66.61 


269 


31 


66.90 


5020053.2 


IRC 00 


"U AT 

b;JN 


2006- 


vo- 


UO -L UU 


50:31 


zuo.oo 


UU.Ul 


269 




uu.yu 


ouzuuoo.o 






9nnfi 

ZUUu- 


05- 


03T20 




268.88 


66.61 


zoy 


30 


66.90 


5020053.4 


IRCOO 


b;N 


2006- 


05- 


03T21 


41:17 


268.88 


66.61 


269 


30 


66.90 


5020053.5 


IRCOO 


b;N 


2006- 


05- 


03T23 


21:02 


268.88 


66.61 


269 


30 


66.90 


5020054.1 


IRCOO 


b;N 


2006- 


04- 


29T20 


18:56 


268.88 


66.61 


269 


35 


66.89 


5020054.2 


IRCOO 


b;N 


2006- 


04- 


30T07 


52:33 


268.88 


66.61 


269 


34 


66.89 


5020054.3 


IRCOO 


b;N 


2006- 


04- 


30T09 


31:47 


268.88 


66.61 


269 


34 


66.89 


5020055.1 


IRCOO 


b;N 


2006- 


05- 


06T00 


52:01 


268.88 


66.61 


269 


28 


66.91 


5020055.2* 


IRCOO 


b;N 


2006- 


05- 


06T05 


49:13 


268.88 


66.61 


269 


27 


66.91 



* MIR-L Data in this observation were not used. 



Table 2. Summary of the survey 



band 


area 


Np* 


T- , t 

integ 




5(7 limit^ 


50% limit'l 




(arcmin^) 




(sec) 




(^Jy) 


(MJy) 


N3 


86.0 


102 


4529 


955 


6.0 


24.0 


S7 


70.3 


300 


4908 


298 


31.5 


47.5 


L15 


77.3 


270 


4417 


277 


71.2 


88.1 



* Number of frames co-added. 

^ Total integration time. 

^ Number of source detected. 

^ Detection limit (sky noise limit). 

II Completeness limit. 




Fig. 1. The locations of the survey fields and the 
SUBARU/Suprime-Cam deep z' band image, superposed on 
IRAS 12 micron map. Blue box and red box are 10x10 
arcmin field of view of the NIR/MIR-S and MIR-L chan- 
nels, respectively. The N3 and S7 fields are covered by the 
Suprimo-Cam image, while the L15 field is partially covered 
by the Suprime-Cam image. The axes shown in the figure 
are in ecliptic coordinate. The field of the AKARI large area 
deep survey (LS NEP DEEP) is indicated as a green circle. 



data obtained from the observations of globular clusters 
in the PV phase (Ita et al. 2007). After the reduction of 
individual frames, we co-added all the frames together. 

In co-adding, we subtracted the self sky image from each 
individual frame which was created by median filtering 
with 20x20 pixel kernel, in order to remove the effect of 
stray light from the Earth shine. Alignment of the frames 
was determined by the position of the point sources in 
each frame in the case of the NIR and MIR-S channels. 
In case of the MIR-L channel, the number of the point 
sources which had sufficient signal-to-noise ratio was not 
sufficient. Therefore, in order to align the MIR-L frames, 
we used the alignment of the MIR-S frames, which were si- 
multaneously observed with the MIR-L in a parallel field, 
using the "coaddLusingS" , an optional task in the IRC 
pipeline. All the 10 pointing data were co-added at once 
by the pipeline in case of the NIR and MIR-S channel. The 
task "coaddLusingS" which was used for the MIR-L data 
does not support a rotation of the FOV and currently can- 
not handle multi-pointing data. Therefore, we co-added 
the MIR-L frames in each pointed observation in order to 
improve the signal-to-noise ratio, and then co-added the 
frames referring to positions of the point sources detected 
in each frame. 

Any spurious events, such as cosmic ray events, mov- 
ing debris, optical and electronic ghost images were re- 
moved at the co-addition stage using a 3 sigma clipping 
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technique. We used the "ccdchp" option (in which the de- 
tector noise performance were assumed) rather than the 
default option "sigclip" (in which the noise performance is 
estimated from the data itself). The average value among 
the frames, instead of the median value, for each pixel was 
computed in order to improve the signal-to-noise ratio of 
the final image. 

The number of frames which were successfully co-added 
were 102, 300 and 270 in the N3, S7 and L15 bands, 
respectively (the L15 band data taken with the target 
ID=5020055.2 was not used since the pipeline failed to 
co-add the L15 data, probably because of large stray light 
from the Earth shine) . One IRC pointing consists of long 
exposure frames together with short exposure frames to 
obtain a wider dynamic range. However, we co-added the 
long exposure frames only. 

Astrometry was based on a comparison of the positions 
of the point sources detected in the AKARI image and 
sources listed in the 2MASS point source catalog. For 
the N3 and S7 band images, a build-in function of the 
pipeline successfully identified the AKARI sources with 
the sources in 2MASS catalog. For the L15 band image, 
automatic matching failed due to the intrinsic red K - L15 
color of the L15 sources. We identified the 2MASS sources 
by eye for the L15 band images. Cross identification with 
the Subaru/Suprimc-Cam deep optical image showed that 
the positional accuracy was better than 2 pixel. 

Finally, we trimmed the area close to the edge of each 
co-added image in which the signal-to-noise ratio was 
poor. The area of the final images are 85.95, 70.27 and 
77.29 arcmin^, for the N3, S7, and L15 bands, respec- 
tively. Figures 2, 3 and 4 show the final images in the N3, 
S7, and L15 bands, respectively. 

4. Source extraction and photometry 

Sextractor (Bertin & Arnouts 1996) was used for source 
extraction. The extraction criterion was a connection of 
five pixels having more than a 1.6532 sigma significance 
above local sky fluctuation. If the flux is equally dis- 
tributed over these five pixels, it corresponds to a 5 sigma 
detection. The actual PSF has a weak concentration at 
its center, and the above criterion gives a slightly higher 
significance than a 5 sigma detection for a point source. 

The flux of each extracted source was evaluated 
by aperture photometry calculated in Sextractor. 
Elliptical aperture photometry with variable size 
(Sextractor's MAGAUTO) was used for the further 
analysis. Parameters of "Kron factor" and "minimum 
radius" were set to the default values of 2.5 and 3.5, 
respectively. The actual apertures used for each detected 
source are overlaid on the AKARI/IRC images in the 
flgures 5, 6 and 7. Magnitude zero point was derived 
by observations of standard stars, and is described in 
Lorcnte ct al. (2007). 

Simple aperture photometry (Sextractor's MAGAPER) 
was also calculated in order to estimate the amount of 
aperture correction. The aperture radius was set to the 
same as used in the standard star observations, 10 and 




Fig. 5. Close up view of the AKARI/IRC N3 band imago. 
Contours are 1, 3 and 10 sigma levels of the local background 
fluctuation. The green circles indicate the positions of de- 
tected sources. The MAGAUTO apertures (cyan) are super- 
posed on each AKARI source. The blue circles indicate the 
apertures on blank sky positions which were used for the cal- 
culation of the detection limit. 




Fig. 6. Close up view of the AKARI/IRC S7 band imago. 
Contours arc 1, 3 and 10 sigma levels of the local background 
fluctuation. The green circles indicate the positions of de- 
tected sources. The MAGAUTO apertures (cyan) arc super- 
posed on each AKARI source. The blue circles indicate the 
apertures on blank sky positions which were used for the cal- 
culation of the detection limit. Some of peaks which have not 
detected as a source by our simple detection criterion have 
counter parts in the deep optical image. Further optimiza- 
tion of detection criterion may improve the detection limit. 




Fig. 2. The final N3 co-added image is shown in equatorial coordinates (rotated 45 degrees counter-clockwise). The green box 
indicates the region over which point source extraction was performed. The size of region is 9.19 x 9.19 arcmin. 




Fig. 3. The final S7 co-added image is shown in equatorial coordinates (rotated 45 degrees counter-clockwise). The green box 
indicates the region over which point source extraction was performed. The size of region is 7.76 x 8.96 arcmin. 




Fig. 4. The final L15 co-added image is sliown in equatorial coordinates (rotated 45 degrees counter-clockwise). The green box 
indicates the region over which point source extraction was performed. The size of region is 8.77 x 9.21 arcmin. 
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Fig. 7. Close up view of the AKARI/IRC L15 band im- 
age. Contours are 1, 3 and 10 sigma levels of the local back- 
ground fluctuation. The green circles indicate the positions 
of detected sources. The MAGAUTO apertures (cyan) are 
superposed on each AKARI source. The blue circles indicate 
the apertures on blank sky positions which were used for the 
calculation of the detection limit. 

7.5 pixels in the NIR and MIR-S/L images, respectively. 
Figure 8 shows the comparison between the flux estima- 
tions by MAGAUTO and MAGAPER for each source. We 
did not apply any aperture correction for further analysis 
because the results using MAGAUTO and MAGAPER 
were reasonably consistent. 

5. Estimation of sensitivity and completeness 

5.1. detection limit (sky noise limit) 

The sensitivity of the point source extraction was esti- 
mated by measuring the fluctuation of the photometry on 
blank sky at random position (sky noise). The positions 
which were close to a source were not used for the mea- 
surement in order to avoid source contamination effects. 

We made simple aperture photometry (a no weight- 
ing, no PSF fitting technique was applied, and no cen- 
tering was performed) at each random position using the 
IRAF/PHOT package. The size of aperture radius was 
set to 1.5 pixel for all three bands. The aperture had 
slightly wider area (7.0 pixels) than the minimum number 
of connected pixels in the extraction criterion (5 pixels), 
a slightly worse estimation of the sensitivity resulted from 
this choice of aperture size. The actual apertures which 
were used for above measurement arc shown in the fig- 
ures 5, 6 and 7. 

In order to evaluate the fluctuation of photometry, the 
histogram of the photometry was fit with a Gaussian func- 
tion. We define the one sigma level of the fluctuation 
as the standard deviation of the Gaussian. Finally, tak- 



ing into account the aperture correction, we calculate a 5 
sigma detection limit for a point source. 

The 5 sigma detection limits arc estimated to be 11.2, 
71.3 and 201 /iJy for the co-added images obtained in 
a single pointed observation, and 6.0, 31.5 and 71.2 fj.3y 
for that in the ten pointed observations, in the N3, S7, 
L15 bands, respectively. The results of the detection limit 
estimation are summarized in Table 4. 

The improvements from one to ten pointings were a 
factor of 1.87, 2.27 and 2.82 in the N3, S7, and L15 bands, 
respectively. These factors in the N3 and S7 bands are 
slightly worse than what is expected from the increase of 
number of observations (i.e., square root of ten). 

5.2. completeness correction 

The completeness of the source extraction was esti- 
mated via Monte Carlo simulations. We added artifi- 
cial sources from a truth catalog into the final co-added 
images, and then extracted the sources again using the 
same extraction parameters. We, then, compared the po- 
sition and magnitude of input sources with the extracted 
sources. If there was a extracted source within 2 pixels 
around an input source and if the difference in magnitude 
between the input source and the extracted source was less 
than 0.5 magnitude, we regarded the event as a successful 
extraction of an artificial truth source, and counted it. 

We define the completeness of the source extraction as 
the number of successfully extracted truth sources divided 
by the total number of artificial input sources that were 
originally added into the original image. 

The calculation was performed by dividing our flux 
range into magnitude bins of size 0.2 magnitudes (i.e. uni- 
form bin size in log Jy) and populating each bin with 20 ar- 
tificial sources. The magnitude distribution of the sources 
was normalized to a flat Euclidian universe in each bin. 
These 20 sources were then added in random positions in 
the image, assuming a circular PSF. In order to avoid self- 
confusion effects among the artificial sources, we did not 
insert a source if the minimum distance to another artifi- 
cial source was less than 20 pixels. Note that the number 
of random sources that can be input at any one time is 
constrained by the array size of the IRC (NIR: 412x512 
pixels, MIR-S/L:256x256 pixels). 

In order to improve the statistical significance of the 
Monte Carlo simulations, we performed 20 statistically 
independent calculations for each magnitude bin, creating 
an effective total of 400 artificial sources. The simulation 
was then subsequently carried out for each magnitude bin 
in turn. 

In order to estimate the error in the calculation of the 
completeness, we then performed 25 statistically indepen- 
dent sets of the above simulations. The results are shown 
in the figure 9, and the completeness limits are summa- 
rized in table 4. 
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Fig. 8. Magnitudes of individual stars given by MAGAUTO are compared with that by MAGAPER with the same aperture radius 
as the observations of the standard stars. No significant systematic difference between MAGAUTO and MAGAPER is seen above 
the 5 (T detection limit, although there is large systematic difference at fainter level. 
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Fig. 9. Differential source counts and completeness estimation is shown for the N3, S7, and L15 bands, respectively. The results 
from both a single pointed and ten pointed observations are plotted. 
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Table 3. Source counts and completeness 













N3 








S7 








L15 






mag^B 


logf Jv) 


N* 




CF 




N 


C 


(7 




N 




(J 


(C) 




<15 


8 


> -2.76 


10 








3 

















15 


8-16 





-2.80 


2 


0.986 





007 


1 


0.985 





006 





0.995 





002 


16 


0-16 


2 


-2.88 


5 


0.983 





005 


1 


0.984 





006 





0.992 





005 


16 


2-16 


4 


-2.96 


2 


0.981 





005 


1 


0.980 





005 


1 


0.988 





005 


16 


4-16 


6 


-3.04 


4 


0.976 





006 


1 


0.982 





005 


3 


0.989 





006 


16 


6-16 


8 


-3.12 


7 


0.975 





008 


1 


0.979 





008 


2 


0.987 





006 


16 


8-17 





-3.20 


3 


0.966 





009 


2 


0.975 





007 


5 


0.983 





005 


17 


0-17 


2 


-3.28 


9 


0.961 





009 


1 


0.972 





009 





0.973 





005 


17 


2-17 


4 


-3.36 


13 


0.946 





009 


3 


0.969 





007 


2 


0.958 





006 


17 


4-17 


6 


-3.44 


10 


0.943 





010 


3 


0.966 





008 


9 


0.940 





009 


17 


6-17 


8 


-3.52 


12 


0.932 





010 


2 


0.959 





009 


8 


0.911 





014 


17 


8-18 





-3.60 


17 


0.929 





010 


3 


0.943 





Oil 


14 


0.874 





018 


18 


0-18 


2 


-3.68 


19 


0.907 





015 


7 


0.927 





012 


11 


0.822 





020 


18 


2-18 


4 


-3.76 


15 


0.888 





015 


7 


0.898 





015 


20 


0.770 





023 


18 


4-18 


6 


-3.84 


25 


0.878 





016 


12 


0.874 





016 


24 


0.692 





016 


18 


6-18 


8 


-3.92 


21 


0.854 





019 


18 


0.822 





013 


19 


0.627 





022 


18 


8-19 





-4.00 


34 


0.829 





019 


19 


0.767 





022 


28 


0.554 





023 


19 


0-19 


2 


-4.08 


43 


0.801 





021 


19 


0.711 





023 


21 


0.480 





026 


19 


2-19 


4 


-4.16 


34 


0.775 





021 


16 


0.635 





026 


35 


0.401 





016 


19 


4-19 


6 


-4.24 


35 


0.734 





024 


24 


0.566 





028 


24 


0.307 





025 


19 


6-19 


8 


-4.32 


37 


0.696 





021 


30 


0.505 





029 


13 


0.217 





022 


19 


8-20 





-4.40 


44 


0.657 





016 


24 


0.432 





021 


17 


0.139 





014 


20 


0-20 


2 


-4.48 


53 


0.601 





021 


21 


0.328 





022 


11 


0.078 





018 


20 


2-20 


4 


-4.56 


46 


0.552 





030 


37 


0.225 





018 


5 


0.037 





010 


20 


4-20 


6 


-4.64 


49 


0.485 





027 


19 


0.145 





020 


4 


0.019 





006 


20 


6-20 


8 


-4.72 


53 


0.426 





028 


12 


0.086 





014 





0.007 





005 


20 


8-21 





-4.80 


65 


0.365 





031 


3 


0.046 





010 


1 


0.003 





003 


21 


0-21 


2 


-4.88 


53 


0.303 





024 


4 


0.020 





006 




0.001 





001 


21 


2-21 


4 


-4.96 


55 


0.243 





019 


2 


0.010 





005 




0.000 





001 


21 


4-21 


6 


-5.04 


48 


0.165 





020 


2 


0.003 





002 




0.000 





000 


21 


6-21 


8 


-5.12 


46 


0.094 





012 




0.001 





002 




0.000 





000 


21 


8-22 





-5.20 


25 


0.052 





014 




0.000 





001 




0.000 





000 


22 


0-22 


2 


-5.28 


26 


0.029 





009 




0.000 





001 




0.000 





000 


22 


2-22 


4 


-5.36 


14 


0.014 





006 




0.000 





000 




0.000 





001 


22 


4-22 


6 


-5.44 


10 


0.008 





004 




0.000 





001 




0.000 





001 


22 


6-22 


8 


-5.52 


7 


0.004 





003 




0.000 





000 




0.000 





001 


22 


8-23 





-5.60 


3 


0.002 





002 




0.000 





000 




0.000 





000 


23 


0-23 


2 


-5.68 


1 


0.001 





001 




0.000 





000 




0.000 





000 



* Number of source detected in each magnitude bin. 
Average of completeness estimation in 25 Monte-Carlo simulations. 
Standard deviation of completeness estimation for the 25 Monte-Carlo simulations. 



Table 4. Detection limit (sky noise limit) and Completeness limit in /ijy 



filter 


N3 


S7 


L15 


limit 


completeness 


detection 


completeness 


detection 


completeness 


detection 


criterion 


30% 


50% 


80% 


ha 


30% 


50% 


80% 


5cr 


30% 


50% 


80% 


5cr 


Ipt 


16.3 


27.7 


81.9 


11.2 


50.8 


75.5 


171.4 


71.3 


93.6 


134.2 


289.6 


200.9 


lOpt 


13.1 


24.0 


83.0 


6.0 


31.7 


47.5 


112.5 


31.5 


57.0 


88.1 


195.2 


71.2 


Ipt/lOpt 


1.24 


1.15 


0.99 


1.87 


1.60 


1.59 


1.52 


2.27 


1.64 


1.52 


1.48 


2.82 
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6. discussion 

6.1. improvement of sensitivity from Ipt to lOpt 

Simple statistics expect a factor of y/lO/1 improvement 
in sensitivity for 10 pointed observations compared to that 
for a single pointed observation. However, table 4 shows 
that the sensitivity did not improve as well as expected. 

One of the effects which degrades the sensitivity at the 
completeness limit is source confusion. In order to esti- 
mate the effect of source confusion, we calculate the source 
number density in the unit of source per beam(Table 5). 
In the case of the N3 band, the density is nearly equal 
to 1/30, the classical limit of source confusion (Condon 
1974; Hogg 2001). Clearly, the observations are limited 
by the source confusion, and this is the reason why there 
is almost no improvement in the sensitivity from 1 to 10 
pointed observations. In the case of the S7 and L15 bands, 
the density is still lower than the classical confusion limit, 
and there is still some room for improvement of the sen- 
sitivity by further pointed observations. 

The detection limit in the L15 band shows a factor of 
3 improvement which is almost as expected from simple 
statistics, while the detection limits in the N3 and S7 show 
somewhat smaller improvements than expected. The de- 
tection limit is estimated by the fluctuation of the local 
background where there are no sources detected. So, a 
smaller improvement in the detection limit means that 
the signal to noise ratio of the sky level is not improved 
by multiple observations, or alternatively that there are 
real fluctuations in the sky. 

There are two possible causes which suppress the 
improvement in the S/N ratio of the sky level. (1) 
Incomplete sky subtraction in the co-addition stage in the 
pipeline creates additional noise in the final image. There 
are two problems which may affect the sky subtraction. 
The first is stray light from the Earth shine, which was 
not uniform over the FOV. The other is scattered light 
from the edge of the detector. Both create non-uniformity 
in the background and may prevent improvement by co- 
adding the frames. In order to eliminate these effects, 
self sky subtraction was performed, i.e. we subtracted the 
median value of surrounding 20x20 pixels from the value 
of each pixel, but there may still be some residual. The 
modeling of the stray light from the Earth shine is still on- 
going and any data taken in the northern sky during the 
months of May, June and July, may suffer from the Earth 
shine problem, and must be handled with great care. We 
have modeled the scattered light from the edge of the de- 
tector, and this model will be included in the next release 
of the IRC imaging "pipeline". (2) Incomplete dithering 
in position in multiple observations. In order to elimi- 
nate systematic errors such as flat fielding error, redun- 
dant observations, in which observations are performed by 
independent detector pixels, are necessary. For this pur- 
pose, dithering of the target position was performed. The 
AOT IRCOO which we used in this observation did not 
perform any dithering during each single pointed obser- 
vation. Dithering of the target position must be realized 
among several pointed observations. Since these observa- 



tions were executed in the PV phase, dithering operations 
among pointed observations was not implemented, at that 
time"^. As described in section 2, only a circular change of 
the target position has been performed in the N3 and S7 
band observations, and one directional shifting of target 
position has been performed in the L15 band observations. 
Systematic errors in the co-addition of frames may result 
from this incompleteness in the dithering. We have since 
implemented dithering among later pointed observations, 
and the observations taken after September 2006 do not 
suffer problems in dithering. 

There is some possibility that we have detected real 
sky fluctuations. Fluctuation analysis of the background, 
however, is still underway, and further deeper observations 
(with a smaller PSF) shall be made in order to produce a 
clear answer to this issue. One thing we have to take care 
is that there may be some sources which have significant 
signal to noise ratio, but could not be extracted by our 
detection criterion (see figure 6). Comparison between the 
deep optical image and the AKARI S7 image shows that 
there are optical counterpart for some faint (below the 
detection limit) AKARI sources. The detection criterion 
which is used in this work is very simple and there is 
much room for improvement. Further optimization of the 
detection criterion may improve the detection limit. 

6.2. 3 and 7 fim source counts 

Figures 10 and 11 show the 3 and 7 /im number count 
results of the AKARI/IRC, ISO/ISOCAM and Spitzer 
IRAC surveys. The counts corrected for in-completeness 
are consistent with each other considering the cosmic vari- 
ance. 

The surveys by IRAC reached deeper than those of the 
IRC because of the larger aperture size of telescope and 
smaller size of the PSF. The IRC, however, has a larger 
FOV (4 times larger than that of IRAC) and a larger 
number of filters (6 filters from 2 to 11 /^m), hence has a 
large merit in wide field and multi-band surveys. 

6.3. 15 fim source counts 

Figure 12 shows the 15 number count results of 
the AKARI/IRC, ISO/ISOCAM and SpitzerlRS surveys. 
This work is currently the widest survey below 100 /iJy 
at this wavelength. The IRC, ISOCAM and IRS resuhs 
are consistent with each other for sources brighter than 
0.2 mJy, and the IRC results confirm the existence of the 
excess compared to the no-evolution source count mod- 
els around 0.4mJy found in the ISO surveys (Elbaz et 
al. 1999). This excess makes itself apparent at around 
the ImJy level. However, there is still a large unccr- 

The specification of the absolute pointing accuracy of the 
AKARI telescope is 30 arcsec. We therefore expect that dither- 
ing in position between multiple pointed observations is achieved 
by scattering from positional error. Actual absolute pointing 
accuracy, in repeated observation for a single target position, 
is measured by this observation, and is turned out to be bet- 
ter than that of expected (the error is within FWHM of the 
PSF). After this observations, we implemented dithering for 
multi-pointed observation by shifting the target position for each 
pointing. 
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Table 5. source number density in the unit of source per beam 



filter 


the number 


of source per beam* 


FWHM of PSF 




1 pt 


lOpt 


in pixel in arcscc 


N3 


1/36.6 


1/31.8 


2.9 4.2 


S7 


1/76.7 


1/56.3 


2.2 5.2 


L15 


1/81.7 


1/61.0 


2.3 5.5 



* The radius of the beam is set to FWHM/2.35 of the PSF. 
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Fig. 10. AKARI/IRC 3 ^tm differential source counts (both raw and corrected for in-completeness) are plotted with the surveys 
carried out by Spitzer/IRAC results. The source counts are normalized to a Euclidian flat universe. Only flux density level above 
50% completeness limit of our results is shown. IRAC channel 1 (3.6 ^m) results (corrected for in-completeness) are from Fazio et 
al. (2004). 
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Fig. 11. AKARI/IRC 7 lira differential source counts (both raw and corrected for in-completeness) are plotted with the surveys 
carried out by ISO/ISOCAM and Spitzer/IRAC results. The source counts are normalized to a Euclidian flat universe. Only flux 
density level above 50% completeness limit of our results is shown. IRAC channel 3 (5.6 lira) results (corrected for in-completeness) 
are from Fazio et al. (2004). ISOCAM LW2 (6.7 iim) results (raw data) are from Sato et al. (2003). 



tainty around the 1-lOmJy flux range due to cosmologi- 
cal variance and cahbration issues with the ISOCAM data 
(Vaisanen et aL 2002). The resuhs from the deep and wide 
field near- and mid- infrared surveys by AKARI/IRC (the 
LS NEP DEEP and LS NEP WIDE surveys; Matsuhara 
et al. 2006) will cover this flux range and will be able to 
settle the issue of the normalization of the bright source 
counts. 

The large discrepancy between the actual source counts 
and predictions by no-evolution models cannot be ex- 
plained singly by K correction effects caused by the strong 
PAH emission from redshifted galaxies and some new pop- 
ulation is required to explain the excess (e.g. Pearson 
2001, Xu et al. 2003, Lagache et al. 2003, Pozzi et al. 
2004, and Pearson 2005). Moreover, evolutionary models 
which fit both the ISO 15 /im and Spitzer 24 y^tm counts 
simultaneously have been already presented by Pearson 
(2005). 

Figure 13 shows the faint end of 15 /im source counts. 
The raw data of the IRC results at the faint end (below 
100/iJy), are consistent with the results from the Spitzer 
IRS le^im PUI surveys in the GOODS fields (Teplitz et al. 
2005a; Teplitz et al. 2005b) in which no completeness cor- 
rection is applied. The source counts corrected for incom- 
pleteness, however, agree better with the faint results from 



the ISO/ISOCAM lensing survey by Altieri et al. (1999), 
and show a slightly higher number density than the pre- 
diction from the model fit to the the ISOCAM results 
(Pearson 2005). The faint end of the source counts, how- 
ever, may suffer from systematic bias such as Eddington 
bias (Hogg & Turner 1998). An ultra deep survey by the 
AKARI/IRC or the Spitzer/IRS W^m PUI is essential to 
produce a definitive answer to the faint source counts. 

7. Summary 

The first near- and mid- infrared deep survey with the 
Infrared Camera onboard AKARI has been successfully 
completed in the performance verification phase of the 
mission. A total of 10 pointed observations have covered 
86.0, 70.3, and 77.3 arcmin^ area, with effective integra- 
tion times of 4529, 4908, and 4417 seconds at 3, 7, and 
15 /xm, have detected 955, 298 and 277 sources, respec- 
tively. The 5 a detection limits of the survey are 6.0, 31.5 
and 71.2/iJy and the 50 % completeness limit are 24.0, 
47.5, and 88.1/xJy at 3, 7, and 15 /im, respectively. The 
observation is limited by source confusion at 3 /im. 

We have confirmed the turnover of the differential 
source counts around 400//Jy discovered by the ISO 15 /im 
surveys. The faint end of our 15 /xm raw source counts 
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Fig. 12. AKARI/IRC 15 fim differential source counts (both raw and corrected for in-completeness) are plotted with the surveys 
carried out by IRAS (shifted from 12 fira), ISO/ISOCAM and Spitzer/IRS results. The source counts are normalized to a Euclidian 
flat universe. A no-evolution model and the contemporary evolutionary model of Pearson (2005) are also plotted for reference. Only 
flux density level above 50% completeness limit of our results is shown. Spitzer/IRS results are from GOODS-N (Teplitz et al. 2005a) 
and GOODS-S (Teplitz et al. 2005b) survey. Data points ISOCAM-DEEP include the ISO cluster lens survey (Altieri et al. 1999), 
and the ISO Hubble Deep Field surveys (Oliver et al. 1997; Aussel et al. 1999). ELAIS-S and ELAIS-N are from the European 
Large- Area ISO Survey (Gruppioni et al. 2002; Serjeant et al. 2000). IRAS counts are from Rush, Malkan & Spinoglio (1993). 



are consistent with the Spitzer IRS results of Teplitz et 
al. (2005a) and the predictions from contemporary evolu- 
tionary models. 

This work is based on the observations in the perfor- 
mance verification phase of the AKARI satellite, which 
is developed and operated by ISAS/JAXA with collab- 
oration with ESA, Universities and companies in Japan, 
Korea, UK, and The Netherlands. The authors would like 
to thank all those who have been involved the AKARI 
project for many years. 
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